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TESTING THE ORIGIN OF HIGH-ENERGY COSMIC RAYS 

A. E. VLADIMIROv', G. JOHANNESSON^, I. V. MOSKALENKO'-^, AND T. A. PORTEr' 

ABSTRACT 

Recent accurate measurements of cosmic-ray (CR) protons and nuclei by ATIC-2, CREAM, and R\MELA 
reveal: a) unexpected spectral hardening in the spectra of CR species above a few hundred GeV per nucleon, 
b) a harder spectrum of He compared to protons, and c) softening of the CR spectra just below the break 
energy. These newly-discovered features may offer a clue to the origin of the observed high-energy Galactic 
CRs. We discuss possible interpretations of these spectral features and make predictions for the secondary CR 
fluxes and secondary to primary ratios, anisotropy of CRs, and diffuse Galactic 7-ray emission in different 
phenomenological scenarios. Our predictions can be tested by currently running or near-future high-energy 
astrophysics experiments. 

Subject headings: astroparticle physics — diffusion — elementary particles — cosmic rays — ISM: general — 
dark matter — diffuse radiation — gamma rays: ISM — infrared: ISM — radio continuum: 
ISM — X-rays: ISM 



1. INTRODUCTION 

The spectrum of cosmic rays (CRs) has offered few clues 
to its origin so far. The only features observed are at 
very high an d ultrahigh energies (see, e.g.. Figure 1 in 

fordvll2 00ll): the so-c alled knee at a few times 10^'^ eV 
ilikov & Khristiansenlll958t iHaungs et al]l2003h . the sec- 
ond "knee" at ~10^^ eV, the "ankle" at higher energies 
(|Abbasi et al. 2005), an d a spectral steepen ing above 10^" 
eV ( Abbas i et al.i i2009l [Abraham et al.l 12010 ). Because of 
the limited size of Galactic accelerators and strength of mag- 
netic fields in the acceleration region (e.g., in supernova rem- 
nants [SNRs]), it is believed that the CRs below the knee 
are Galactic, while above the knee they have an extragalac- 
tic origin, with the knee itself being due to propagation ef- 
fects and a transition between the two po pulations of CRs 
(iBerezinskii et al.l ll990: Stron g et al.ll2007l) . 

The power-law spectrum below the knee is thought to 
be the result of C R acceleration in SNR shocks (see, e.g., 
iDrury et aH 120011) . which is steepened to the observed in- 
dex ^ 2.75 by propagation in the interstellar medium 
(ISM) and eventual leakage from the Galaxy. The inter- 
stellar diffusion coefficient is typically assumed to be a 
power law in particle rigidity, based on numerous studies 
of m agnetohydrodynamical turbulence (see, e.g., iBiskamp 
120031) . The turbulent cascade often leads to a distribu- 
tion of magnetic energy that is well described by a power 
law. For energies below ^ 20 GeV nucleon^ ^, the CR 
spectrum flattens due to the modulation in the heliosphere 
— a combined effect of the solar wind and heliospheric 
magnetic field. Measurements of CR composition below 
a few GeV nucleon"^ offer detailed informa tion on el- 
emental and isotopic abundances (Engelmann et al.l 119901: 
IWiedenbec k et aLliooTt lObermeie r et al.l2011 ), including the 
peaked shape of the secondary-to-primary nuclei ratio (e.g., 
B/C, sub-Fe/Fe) and abundances of long-lived radionuclides 
(such as i"Be, ^^Al, ^^Cl, and ^"^Mn). These measure- 
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ments are used to derive the (model-dependent) diffusion 
coefficient and the size of the Galactic volume filled with 
CRs dS trong & Moskalenko 1998; Ptus kin & Sou toul 199^ 
IWebber & SoutouLJ998,) . the so-called halo. Models of CR 
prop agation are in reasonable agreement with available data 
(e.g. .'Stron g et al.ll2007l:lTrotta et alj201 Ih . with a few excep- 
tions, including the un expected rise in the positron fraction 
observed by PAMELA (lAdriani et al.ll2009h . 

The data recently collected by three experiments, ATIC- 
2 (Wefel et al. 2008; Panov et al. 2001), CREAM (Ahn et al] 
aOlOEiYoon et al.n2011.) . and PAMELA ( lAdriani et al..i20lC 
indicate a break (hardening) of the spectra of the most abun- 
dant CR species above a rigidity of a few hundred GV. 
The break rigidity, pbr, is best measured by PAMELA and 
occurs at approximately the same rigidity for p and He, 
Pbr = 240 GV. The PAMELA data for 10 GV < p < pbr 
agre e very well with the earlier data from AMS and BESS 
(see Alcaraz et al. 2000; ' Haino etal.1 12004' and Fi gure 1 of 
Adriani et al. 2011), while ATIC-2 data points for p < pbi 
are somewhat lower. We take the PAMELA data as the most 
accurate for p < pbr- For p > pbr, ATIC-2 results agree well 
with those of CREAM. The change in the spectral index (be- 
low/above the break) is estimated as Abr = 7(> Phi ) — j{< 
Phi) — 0.15, and is the same for protons and He. 

Another important feature of the CR spectra discovered by 
these experiments is the difference between the spectral in- 
dices of CR pr otons and He. This h as been speculated for a 
long time (e.g..' Biermann et al]|1995' , and references therein), 
but the experimental uncertainties were too large to be conclu- 
sive (see the collection of CR proton and He measurements 
in Moskalenko et al. 2002). The new measurements by the 
ATIC-2, CREAM, and PAMELA experiments confirm this 
with high significance. The spectrum of He is found to be 
harder than the spectrum of protons for energies up to, at least, 
10"* GeV nucleon^ ^. The difference between the proton and 
He spectral indices calculated by Adriani et al. (2011) using 
the PAMELA data is A7 = 0.10, and it is approximately the 
same above and below pbr- Within the statistical and system- 
atic uncertainty, the measured p/He flux ratio appears to be a 
smooth function of rigidity, continuous at pbr- This shows that 
the difference in the spectral slope of protons and He nuclei 
persists into the ultra-relativistic regime. 
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There is also fine structure in the spectra that may provide 
some clues to the nature of the observed features: PAMELA 
data clearly show a spectral softening at the break rigidity 
(which we refer to as the "dip", below). lAdriani et alJ (1201 Ih 
have shown the softening to be statistically significant at the 
95% confidence level for the spectra as functions of particle 
rigidity, and at the 99.7% level for the same data in terms of 
kinetic energy per nucleon. The softening is more pronounced 
in the He spectrum. 

Rather than proposing a detailed interpretation of the ob- 
served features, in this paper we discuss broad categories of 
models, hereafter called Scenarios, and propose their observa- 
tional tests. A particular reaUzation of each scenario is called 
Calculation. The quantitative analysis is done using the GAL- 
PROP cod^ ("Strong & Moskalenko|[T998l) . 

We study the interpretations of the p/He ratio variation sep- 
arately from the interpretations of the spectral break at pbr- In 
Section im we introduce the reference scenario based on the 
pre-PAMELA data. In Section 12.21 we discuss possible ex- 
planations of the p/He ratio decline with energy: inherent na- 
ture of CR sources and spallation effects. Section l23] presents 
four physical scenarios that could lead to the observed spec- 
tral break at pbi- injection effects, propagation effects, and 
local low- or high-energy CR source. The framewor k of our 
CR propagation calculations is described in Se ctio ns [ 3. II and 
13.21 and specific calculation setups in Section [33] The sec- 
tions following that discuss our results and their implications 
to CR observations and CR propagation modeUng. 

2. SCENARIOS 

2.L Reference Case 

Scenario R: Reference scenario. First, we introduce a ref- 
erence case based on the pre-PAMELA data. In this scenario, 
the CR injection spectrum above 10 GV is a single power law 
up to the "knee" in the CR spectrum, with the same spectral 
index for all CR species. The rigidity dependence of the diffu- 
sion coefficient at high energies is also taken as a single power 
law for all energies. The CR source distribution is described 
in Section |3] 

Scenario R provides reasonable agreement with the pre- 
PAMELA data, but it cannot reproduce the spectral features 
evident in the new data discussed in this paper: the difference 
between proton and He spectra, the spectral break, or the dip. 
Below, we describe several broad categories of models that 
encompass viable explanations for these new features. The 
comparison of predictions for these other models for quanti- 
ties other than CR proton. He, and electron spectra with pre- 
dictions of Calculation R qualitatively illustrates the signifi- 
cance of the difference between different scenarios. 

2.2. p/He Ratio: Acceleration and Spallation Hypotheses 

The confirmation of a significant difference between proton 
and He spectral indices poses a challenge for theories of CR 
acceleration and propagation. Whatever the physical cause of 
this difference in spectra may be, it seems to affect heavier nu- 
clei in the same way as it does He (see, e.g.. lAhn et al.ll2009[ 
for spectra of nuclei), giving them a harder spectrum than that 
of protons. 

Diffusive shock acceleration (DSA) predicts the spectrum 
of He trapped in a shock to be harder than that of protons 

^ The project Web site |http://galprop.stanford.edu7] 



due to its lower Z/A ratio, b ut only for non-relativistic ener- 
gies (e.g.. lEllison et al.lll997l) . However, in particles escaping 
from a shock, p/He ratio may decline with energy, if DSA 
is rapid, and the injection of He into the acceleration pro- 
cess varies in a way that enhances He acceleration in young 
shocks. This could happe n due to the inherent property of par- 
ticle injectio n in shocks dMaUcov et al.ll20"T2l) . or if the abun- 
dance of He dOhira & Iokal201 ll) or magnetic field orientation 
jBiermann et al.ll995b is inhomogeneous in the SNR environ- 
ment. We encompass these mechanisms into the Acceleration 
Hypothesis {Hypothesis A), which is discussed below. Note 
that propagation effects may contribute to the p/He spectral 
difference because the second-order Fermi process (reaccel- 
eration) in the interstellar medium makes the He spectrum 
harder due to its lower Z/A ratio. This effect, howev er, does 
not ex tend to the ultra-relativistic regime (e.g.. Strong et alJ 
l2007h . 

An alternative idea, suggested by iBlasi & Amatol (l2012ah . 
is that that spallation of CR nuclei {Z > 1) may lead to hard- 
ening of their spectra. This is because the lower energy CRs 
have longer confinement times in the Galaxy, and their flux is 
depleted by spallation more than the flux of higher energy nu- 
clei. Hardening occurs only if the spallation timescale is short 
compa red to the confineme nt timescale of the nuclei. Note 
that Blasi & Amatol ( l2012al) consider CRs at energies above 
1 TeV and do not attempt to make their model consistent with 
the CR data at low energies, where various effects, such as 
stochastic reacceleration and significant production of secon- 
daries, come into play. We investigate this idea, extending it 
to lower energies, and hereafter refer to it as Hypothesis S. 
Our calculations include spallation of all nuclei species at all 
energies, the default with GALPROP. However, as our results 
for Calculation R show (see Section HI, the effect of spalla- 
tion on the He spectrum is insignificant, and the p/He ratio 
above 10 GV is flat. Below, we demonstrate that, with some 
model tuning (i.e.. Calculation Si,2), fragmentation may in- 
deed lead to hardening of the He spectrum. We also assess 
the consequences of the required model modifications. 

Note that in this section we refer to the results of our calcu- 
lations for Hypothesis S. The framework and details of these 
calculations are formally introduced later in the text in Sec- 
tion |3] However, we would like to briefly discuss Hypothe- 
sis S here, because in the rest of the paper we discard Hypoth- 
esis S and adopt the assumption that the p/He decline is caused 
by the nature of CR accelerators (Hypothesis A). The reader 
interested in our reasoning for discarding the spallation ef- 
fects hypothesis may find the explanation in this section. And 
a detailed description of our physical model, computational 
method and data sources can be found in Sections |3] and ID 
where we discuss the spectral break at pbi ■ 

Hypothesis S: Spallation effects. The fraction of frag- 
mented CR nuclear species depends on their total inelas- 
tic cross section and the effective grammage encountered by 
the CR species in the Galaxy. Inelastic cr oss section fits 
used i n iBlasi & Amatol (l2012al) . taken from iHorandel et alj 
(120071) , are somewhat larger than those used in our calcu- 
lations (iBar ashenkov 1993; Ba rashenkov & Polvan ski 1994). 
Besides that, the gas number density used in calculations 
by Blasi & Amato (2012a) yields a significantly larger gram- 
mage than in our standard models. To segregate these effects, 
we construct two calculation setups for Hypothesis S: Cal- 
culation Si and Calcu lation S2. Both calcul ations adopt the 
cross section fits from IHorandel et aTl ( 120071) . Calculation S 
uses parameters similar to that of the reference Calculation R, 
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but a sli ghtly smaller diffusion c oefficient to match the value 
used by iBlasi & Amatol (l2012al) . In Calculation S2, we ad- 
ditionally increase the gas number density relative to Calcu- 
lation Rhy a factor of two. Note that these calculations use 
the GALPROP code, which was adapted to incorporate the 
above-mentioned inelastic cross sections, whereas the pro- 
duction of fragments (daughter isotopes) is calculated using 
a standard set of cross sections and remains unchanged. 

Our calculations show that with our standard gas distri - 
bution based on astronomical data (Mo skalenko et al.]l2002h . 
Calculation Si, the amount of hardening is insufficient to pro- 
vide agreement with the PAMELA p/He ratio. The required 
hardening of the He (^He + *He) spectrum can be achieved 
only if we assume a considerable increase in grammage and 
simultaneously adopt a set of total inelastic cross sections 
from Horandel et al. (12007b . Calculation 82- However, this 
leads to an overproduction of secondary species in CRs, such 
as antiprotons and boron, so that the calculated B/C ratio does 
not agree with the data. The p/He ratio obtained in Calcu- 
lation Si and Calculation S2 is shown in the top panels of 
Figure [7] (see Section 13.21 for an overview of the plain diffu- 
sion and diffusive-reacceleration models). The calculated B/C 
ratio and p flux are shown in the middle and bottom panels, 
respectively. Read on and see Figures [lj{3] for details on the 
parameters of the calculations shown in Figure [T] 

Another important point to consider here is that the mea- 
surements of PAMELA, ATIC-2, and CREAM are not sensi- 
tive to the isotopic composition of CR fluxes. The He fluxes 
reported by these experiments and used throughout this paper 
are, in fact, the sum of '^He and '*He species. The dominant 
channel of *He spallation is the reaction 

^Hc + p ^ ^Hc + X. 

This reaction leads to a hardening of the interstellar ^He spec- 
trum because lower energy nuclei experience more spallation 
events. However, due to production of ■^He in the same reac- 
tion, the total He spectrum does not harden as much as ''He 
alone. Further spallation of secondary '^He, as well as frag- 
mentation of '*He into products other than '^He, eventually 
leads to the total He spectrum hardening. Still, the effect of 
spallation on the total He spectrum is not as str ong as on ''He 
alone. Equation (2.1) in lBlasi & Amatol (120 12al) indicates that 
■^He was not included in their calculations. Therefore, their 
results would be relevant only for the p/^He ratio. 

This is illustrated in Figure|7] where we also plot the ratio of 
p/'^He in Calculation S2 for reference. It can be seen that the 
overall shape of the p/'^He ratio matches the measured p/He 
ratio well; however, a significant fraction of secondary ^He 
changes the shape so that the calculated p/He ratio cannot be 
adjusted to match the data simultaneously at all rigidities (1 
GV - 10 TV). 

Hypothesis A: Acceleration effects. We have concluded that 
the adjustments of propagation model required to reproduce 
the observed p/He ratio in Hypothesis S2 conflict with the 
measurements of secondary CR species. Therefore, in the rest 
of this work we adopt an alternative to Hypothesis S, which we 
call the 'Acceleration hypotheses". Hypothesis A represents 
the idea that the nature of CR accelerators is responsible for a 
harder spectrum of He than p (see references at the beginning 
of this section). To incorporate Hypothesis A into our calcu- 
lations, we use an ad hoc modification for the CR injection 
spectra. That is, our calculations assume that nuclei heavier 
than H are injected into the ISM with a harder spectrum than 



protons. The difference between the spectral indices of pro- 
tons and heavier nuclei is the same for all rigidities and is 
Ap/Hc = 0.07. 

We do not present a separate calculation for the Hypothe- 
sis A in this paper. Instead, we incorporate Hypothesis A into 
calculations that study scenarios P, I, L and H for the break in 
the p an d He spectra (those scenarios are introduced in Sec- 
tion |2]3]l. Figure |5] shows the p/He ratio resulting from this 
modification (see below for parameters of other calculations 
shown in this figure). The figure illustrate that the data re- 
ported by PAMELA, ATIC-2, and CREAM can be reproduced 
by including Ap/Ho = 0.07. 

2.3. Spectral Break and Dip: Propagation, Injection and 
Local Source Scenarios 

We consider the following scenarios for an explanation of 
the break at pbi and the dip just below pbi- (1) interstellar 
propagation effects, (2) modification of CR injection spec- 
trum at the sources, (3) composite Galactic CR spectrum, (4) 
effects of local sources at low energies (p < pbr), and (5) 
effects of local sources at high energies {p > pbi )- Particu- 
lar realizations of these scenarios (calculations) are discussed 
in detail in Section 13.31 their parameters are summarized in 
Tables[I]and|2l 

Scenario P: interstellar Propagation effects. Transport 
of CRs in the ISM is subject to considerable uncertain- 
ties, because the properties of interstellar magnetic turbu- 
lence are not very well known (jElmegreen & Scalo 2004|; 
Scalo & Eknegreenl I2004Y This makes CR observations a 
valuable indirect probe of quantitative features of particle 
transport (e.g., the diffusion coefficient, D) in the Galaxy. 
Therefore, in this scenario, the break in the observed pro- 
ton and He spectra is attributed to a change in CR transport 
properties at rigidity pbr- This scenario is represented by Cal- 
culation P, which has a break in the rigidity dependence of 
the diffusion coefficient at p — pbr- For p < pbr, we use the 
functional form of D(p ) obtained in t he ear lier comprehensive 
analysis of CR data bv lTrotta et alj (1201 Ih . and for p > pbr, 
we adjust the rigidity dependence of D{p) to match the ob- 
servations of PAMELA, ATIC-2, and CREAM, as discussed 
above. 

Scenario I (a): CR Injection effects, source with a spec- 
tral break interpreta tion. Existing models of CR production 
by SNR shocks (e.g.JCap rioli et al.ll2010l:IPfiiskin et al.l20Tol) 
predict a smooth spectrum of CR particles injected by an SNR 
into the Galaxy. Such models usually consider a shock in 
a semi-infinite medium or assume spherical symmetry. The 
spectrum predicted by these models may gradually harden 
with energy between 10 GeV and 100 TeV, but not as rapidly 
as in the PAMELA data. Note that particle transport, magnetic 
turbulence generation, and nonlinear feedback of particles and 
magnetic fields on shock structure are not strictly constrained 
in these models. The spectrum of particles leaking from an 
SNR shock has never been observed directly. It is therefore 
conceivable that with some parameter tuning, present mod- 
els of particle acceleration may predict a more pronounced 
hardening in the spectrum of particles injected into the ISM, 
consistent with the new data. Alternatively, particle acceler- 
ation models that take into account the asymmetry of SNRs 
may predict a break in the particle spectru m produced by a 
single SNR. For example, in the model of iBiermann et alj 
(2010), the break, or upturn, occurs due to the contribution of 
the SNR's polar cap. This case, hereafter referred to as Sce- 
nario I (a), is represented by Calculation I, which features a 
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Galaxy- wide source spectrum with a hardening at pbi - The 
diffusion coefficient does not have a break in this scenario. 

Scenario I (b): CR Injection effects, composite source in- 
terpretation. While SNRs (isolated or in superbubble regions) 
are believed to be the primary sources of Galactic CRs, dif- 
ferent classes of supernovae and their environments, as well 
as other CR sources, can combine to produce the observed 
CR spectrum. Generally speaking, different types of CR 
sources could have different spatial distributions throughout 
the Galaxy. In this work, we make the simplifying assumption 
(1) that there are only two types of CR sources, and (2) the 
spatial distributions are the same for both types of CR sources. 
If one source dominates the low energy part of the CR spec- 
trum, and the other the high energy part, then Calculation I, 
with a hardening of the Galactic CR source at pbi , adequately 
encompasses this composite source scenario as well. This 
scenario may be generali zed to a d istribution of CR sources 
with different parameters. I Yuan et a l. (201 1) have shown that 
in general, dispersion in the CR source spectral indices re- 
sults in the concavity of the observed CR spectrum. We use 
the same computational setup to calculate the observed quan- 
tities for Scenario I (a) and Scenario I (6), and we call it just 
Calculation I. A subtle advantage of the composite source in- 
terpretation of Calculation I (i.e., in Scenario I (6)) is its abil- 
ity to explain the dip more naturally than the source with an 
inherent break scenario (see the discussion of the dip in Sec- 
tion |4|. 

Scenario L: local Low-energy source. This scenario encom- 
passes interpretations that assume that the observed spectral 
break is caused by a local source dominating the CR spec- 
trum at low rigidities, p < pbr- Unlike Scenario I (6), the 
present scenario assumes that the low-energy source is not 
typical for the Galaxy as a whole. This scenario is formu- 
lated as Calculation L, in which the Galactic CR spectrum is 
hard, matching the observations of PAMELA, ATIC-2, and 
CREAM for p > pbr- For p < pbi , the flux of Galactic CRs 
is lower than the observed flux, and we assume that the dif- 
ference is accounted for by the hypothetical local source. We 
assume the extreme case of a very local low energy source. 
This means that we do not calculate propagation of CRs from 
that source and only the Galactic sources with the hard spec- 
trum are used to calculate the production of secondaries and 
the diffuse Galactic 7-ray emission. This scenario contrasts 
with Scenario I (&), where the sources of low-energy CRs are 
distributed across the Galaxy. The case of intermediate local 
source extent falls in between Scenario L and Scenario I (b). 

Scenario H: local High-energy source. This scenario is 
analogous to Scenario L, but with Galactic sources dominat- 
ing the CR flux for p < pbi , and the spectral break produced 
by a local high-energy source dominating the observed flux 
for p > pbr- The calculation representing this scenario is re- 
ferred to as Calculation H. The assumption of the high-energy 
source being very local is made in this calculation identically 
to how it was done in Calculation L, i.e., the production of 
secondaries and the diffuse Galactic 7-ray emission is deter- 
mined solely by the Galactic CR sources. 

3. CALCULATIONS 

3.L GALPROPCode 

The GALPROP project began in late 1990s 
dStrong & Moskalenkol [1998) and has been in continu- 
ous development since. The code is available from the 
dedicated Web site where a facility for users to run the 



code via online forms in a web browseiQ is also provided 
dVladimirov et al.ll201 lb . 

The GALPROP code solves the CR transport equation for 
a given source distribution and boundary conditions for all 
CR species. This equation includes diffusion, a galactic 
wind (convection), diffusive reacceleration in the ISM, energy 
losses, nuclear fragmentation, radioactive decay, and produc- 
tion of secondary particles and isotopes: 
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where ip — ip{r, p, t) is the CR number density per unit total 
particle momentum (i.e., ip{p)dp — Anp^ f {p)dp in terms of 
phase-space density f{p)), q{v,p) is the source term, Dxx is 
the spatial diffusion coefficient, V is the convection velocity, 
reacceleration is described as diffusion in momentum space 
with diffusion coefficient Dpp, p = dp/dt is the momentum 
loss rate, t/ is the time scale for fragmentation, and is the 
timescale for radioactive decay. The numerical sol ution of the 
transport equation is based on a Crank-Nicholson (iPress et al.l 
[l992) implicit second-order scheme. The spatial boundary 
conditions assume free particle escape, e.g., ijj{Rh, z,p) = 
ip{R, ±Zh,p) = 0, where Rh and Zh are the boundaries for a 
cylindrically symmetric geometry. 
The source function q is 



p) = gpri(r, p) + V gscc(r, p) 



(2) 



where gpii represents the primary CR sources, and the gscc 
term is for the sources of secondary isotopes (i.e., nuclear re- 
actions in the ISM), and p = pc/Ze is the magnetic rigidity 
where p is momentum and Ze is the charge. The distribution 
of primary Galactic CR sources used in this work is based o n 
the supernova distribution from lCase & Bhattacharyal (11998'). 

While GALPROP's numerical scheme can accommodate 
arbitrary energy dependence of CR source function, in this 
work we parameterize the source function of primary nuclei 
as a broken power law in particle rigidity: 



9pri(r,p) oc Q(r) 



(3) 



( go for p < pinjo 

9=1 91 for pinjo < P < Pinjl , 
I 92 for p > Pinjl 

The source function for primary CR leptons is similar to that 
of nuclei, but with up to two breaks. 
Likewise, the spatial diffusion coefficient is given by 

s 

(4) 



( 60 for p < po 
5=1 ^1 for po < /O < Pi , 
[ 82 for p> pi 

where Do is the normalization at rigidity po and /3 = v/c. 
The power-law index S = 1/3 corresponds to Kolmogorov 
diffusion (see Section lsTl i. 

http://galprop.stanford.edu/webrun/ 
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GALPROP solves the time-dependent propagation equa- 
tion, and in this work, the steady-state solutions of Equa- 
tion ^ are obtained assuming that the source functions are 
time-independent and integrating the equation over a long 
enough time interval. The accelerated solution technique is 
used, where the initial time step. At = 10^ yr, is large com- 
pared to the propagation timescale, and after Ng — 20 itera- 
tions. At is reduced by a factor of two, etc., until At becomes 
small compared to the shortest timescale in the system (in our 
case, 10 yr, to accommodate the rapid energy losses of lep- 
tons). 

The details of physical processes and data used in the GAL- 
PROP code, as well as the numerical scheme, can be found 
elsewh ere. A complete l i st of r elevant publications is avail- 
able in IVladimirov et al.l (1201 1\. the aforementioned GAL- 
PROP Web site contains additional information and publica- 
tions. 

3.2. Dijfusive-Reacceleration and Plain Diffusion Models 

Previous studies have shown that the available CR data can 
be explained in one of the two common propagation models: 
the diffusive- reaccele ration (D-R) model and the plain diffu- 
sion (PD) model. These models have been used in a number of 
studies utilizing the GALPROP code (e.g., Moskalenko et a^ 

etall 



20021; IStrong et all [2004; Ptuskin et al.. .2006^ Abdo ( 



20091 and references therein). 



The D-R model assumes that low energy CRs in the ISM 
participate in the second-order Fermi acceleration process. 
This process is believed to be caused by stochastic collisions 
of CR particles with moving magnetic structures. Averaged 
in time, such collisions result in particle diffusion in momen- 
tum space with a diffusion coefficient Dpp, which increases 
the mean energy of low energy parti cles. If reacceleration 
is incl uded, Z)j,i, is related to D^x (iBerezinskii et al.l 119901: 
lSeo&Ptuskia,1994.) : 



DppDxx 



3(5(4 -52)(4-5)u.' 



(5) 



where w characterizes the level of turbulence (we take w = 
1 because only the quantity v\^^/w is relevant), and 5 — 
1/3 for a Kolmogorov spectrum of interstellar turbulence 
(Kolmogorov 1941) or 5 = 1/2 for a Kraichnan cascade 
dlroshnikov 1964: Kraichnanll 1 965h . but can also be arbitrary. 
Matching the B/C ratio below 1 GeV in D-R models is known 
to require large values of WAif ■ In order to avoid a large bump 
in the proton spectrum at low energies, the D-R model re- 
quires a break in the CR injection function around p = 10 GV. 

The PD model assumes no reacceleration process, which 
corresponds to vau = 0. No break in the CR injection func- 
tion is required to fit the proton and helium spectra at low en- 
ergies in the PD model. However, in order to fit the B/C data 
below 1 GeV nucleon^^, the PD model requires a low-energy 
break in the diffusion coefficient. Specifically, the diffusion 
coefficient in the PD model must decrease with increasing 
energy below 4 GV in order to fit the B/C measurements be- 
low 1 GeV nucleon^^. A possibl e physical justificatio n of 
such behavior of D{p) is given bv Ptuskin et al.l (I2006cl) and 
involves turbulence dissipation in the ISM. 

Each of our calculations {Calculation R, P, I, L, and H) is 
presented in two versions: one for the D-R and another for the 
PD model. 

3.3. Calculation Setups 



The parameters of our calculations are summarized in Ta- 
bles[T]and|2l Figures [T}{3] show the diffusion coefficients and 
the injection spectra used for the different scenarios. 

Calculation R is the reference case for this study. The list 
below outlines the key parameters of this calculation for the 
D-R and PD models. 

1. For Calculation R in the D-R model, we chose go — 
— 1.9, gi = g2 = —2.4 (i.e., no break at pinji) and 
Pinjo = 11 GV for all nucleons, which is consis- 
tent with the findings of Trottaet al. (2011) based on 
a Bayesian analysis of a D-R model using the GAL- 
PROP code. For the PD model, a low energy break 
in the injection spectrum is unnecessary, and our Cal- 
culation R uses 5o = 5i = 52 = —2.1. The elec- 
tron injection spectrum for Calculation R in the D-R 
model is similar to that from Ackermann et al.' (201(f), 
with spectral index in rigidity r=l .60/2. 50 below/above 
a break rigidity of Pco=^ GV, and a second steepening 
to r=5.0 above Pc2=2 TV. For the PD model, we use 
r=2.35 above 4 GeV. Note that the low-energy break in 
the electron injection function can be explained without 
the corresponding break in the proton injection spec- 
trum, because electrons and protons may be accelerated 
in sources via different mechanisms. 

2. For the diffusion coefficient. Calculation R uses S — 
0.30 in the D-R model, and S = 0.60 in the PD model. 
The normalization for the diffusion coefficient in the D- 
R model is Do = 5.75 x lO^* cm s^^ at p = 4 GV, 
wh ich is consi s tent w ith the best-fit values obtained 
by iTrotta et al] ( 1201 Ih . For the PD model, we use 
Do = 3.0 X 10^® cm s~\ in order to match the B/C 
observations of HEAO, TRACER, and CREAM (see 
Section |4~4l l. In this work, we construct the PD model 
with a constant, rather than decreasing, diffusion co- 
efficient below 4 GV, even though the model does not 
reproduce the low energy B/C data. This is done in or- 
der to illustrate the possibility that a local low-energy 
CR source (Scenario L) can simultaneously explain the 
break in p. He spectra and fit the low-energy B/C data 
without requiring the diffusion coefficient to decrease 
with increasing energy (see below). 

3. Finally, to model reacceleration in the D-R model, we 
chose VAu = 32 km s~^, the halo size — 4 kpc, and 
the normalization of the propagated CR proton spec- 
trum was tuned to the observed flux Np — 10.7 x 
10^12 sr-i MeV-i for p = 10^ GV. These 
values were obtained by slightly adjust ing the best- 
fit values obtained by 'Trotta etal] (1201 lb to achieve a 
good agreement with the PAMELA proton spectrum for 
P < Pinjo- These adjusted values are still within one 
mean squar e devia tion of the posterior mean found by 
ITrotta et alj dloTlh . In all fi gures, black lines represent 
the input and output quantities pertaining to Calcula- 
tion R. 

Calculation P has the same parameters as the reference 
Calculation R, except: (1) the injection spectrum of protons 
(electrons) above the low energy break pinjo (Pco) has a softer 
power-law index to give agreement with the PAMELA data 
below pbr; (2) the injection spectra of He and heavier ele- 
ments (A > 1) have a power-law index harder than that of 
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protons by Ap/Hc for all rigidities (this represents Hypothe- 
sis A); (3) the rigidity dependence of the diffusion coefficient. 
Equation (|5]l, has a break at rigidity pi (i.e., 5i ^ (^2). The 
break in the rigidity dependence of the diffusion coefficient 
is introduced to match the observed break in the CR spec- 
trum. Besides, we choose pi k, p-^^ so that pi is slightly 
larger than pbr for better agreement with the data for p > pbr- 
The normalization of the proton flux has been adjusted in Cal- 
culation P, along with the abundance of He, to agree with 
PAMELA data at all energies (the abundances of heavier nu- 
clei were changed by the same factor as He). The results for 
Calculation P are shown with blue lines in all figures. 

Calculation I differs from Calculation R in the following 
ways: (1) the index of the proton injection spectrum is softer 
than in Calculation R for pinjo < P < Pinji; (2) the injec- 
tion spectrum has two breaks, i.e., gi ^ gi (this represents 
Scenario I (a) and Scenario I (b)); (3) electrons also have a 
softer spectrum for pco < P < Pci, and a break at pd, and 
(iv) nuclei are injected with a harder spectrum than protons 
(Hypothesis A). This calculation produces a CR spectrum at 
Earth with a break at pbr closely matching that of Calcula- 
tion P, but due to a different physical assumption. Namely, 
it is the spectral break in the CR injection spectrum that pro- 
duces the break in Calculation I, whereas in Calculation P, it 
occurs because of a break in the diffusion coefficient. Note 
that the high-energy break in the electron injection spectrum 
must be stronger than for protons, and it must occur at a lower 
rigidity than in protons, in order to obtain agreement with the 
electron spectra observed by the FermZ-LAT and PAMELA. 
The results of this calculation are shown as green lines. 

Calculation H combines two components. One component 
is produced by the Galactic CR sources and propagated using 
GALPROP with the same parameters as Calculation I. This 
component does not have a break in the CR injection spec- 
trum at pinji and its normalization was tuned to match the 
proton and He spectra for p < pinji- Another component is 
produced by a hypothetical local source, which contributes to 
the total flux only for p > pbi ■ We do not calculate CR propa- 
gation for the local source; instead, we calculate its spectrum 
at Earth by subtracting the Galactic source spectrum from the 
data of ATIC-2 and CREAM for p > pbr- As discussed in 
SectionlH this represents the assumption that the local source 
is nearby and not very powerful. To compute secondary par- 
ticles and isotopic ratios in this calculation, we assume that 
isotopic abundances in the local and Galactic sources are sim- 
ilar, and that the local source supplies no secondary particles 
at Earth. This lowers, for instance, the B/C and p/p ratios for 
p > Phi- The diffuse 7-ray emission from the Galaxy is cal- 
culated using only CR fluxes from the Galactic source. Gray 
lines represent this calculation in plots. 

Calculation L features a local source that contributes to the 
Zow-energy part of the CR spectrum (p < pbr)- The local 
source is included in the same way as in Calculation H, i.e., its 
propagated spectrum is calculated as the difference between 
the observed CR spectrum and the propagated Galactic com- 
ponent. As in Calculation H, we do not calculate the propa- 
gation of CRs from the local source and assume that its flux 
contains no secondary species. However, this calculation is 
very different from all other, because the low-energy CRs are 
a mix of particles which have undergone Galactic propaga- 
tion and recently accelerated particles from the local source. 
Because of that the propagation parameters for this calcula- 
tion should be estimated simultaneously with the parameters 
of the local source, as in .Moskalenko et al.. (.2003i) . 



Indeed, assuming that CRs from the local source are pro- 
duced so recently that they contain no secondary nuclei, and 
that there is no primary boron in CRs, one can calculate 
the B/C ratio for Scenario L in the energy range where the 
local source contribution is non-negligible. Matching the 
B/C ratio in the D-R is quite challenging because the lo- 
cal source contribution reduces the B/C ratio in the range 
1-10 GeV nucleon^^. This reduction may be compen- 
sated by assuming a lower diffusion coefficient in this en- 
ergy range. At the same time, the diffusion coefficient above 
10 GeV nucleon^^ cannot be changed very much in order to 
maintain agreement with high-energy B/C data. The above 
considerations necessitate a larger value of Si. With a greater 
Si, the diffusive reacceleration of low-energy protons be- 
comes too strong, resulting in an overprediction of the proton 
flux around 1 GeV, and thus UAif should be reduced. In fact, 
we found that in this scenario, the best agreement with the 
PAMELA data for protons is achieved with vau — (i.e., no 
reacceleration) and So=0 below a 4 GV (i.e., the PD model). 
This is because any finite vau hardens the proton spectrum 
below 1 GeV too much to match the PAMELA data. 

Considering the above, we chose to keep Calculation L in 
the D-R model with unchanged diffusion coefficient, in or- 
der to illustrate the problem in the B/C fitting. And in the 
PD model, the reduction of the B/C ratio below 1 GeV was 
beneficial for agreement with data, because in all other PD 
calculations, B/C was overpredicted. Quantitatively matching 
the ACE data in the PD model requires the flux of the local 
source to be relatively small. This dictates our choice of a 
concave Galactic source spectrum for the Calculation L in the 
PD model, i.e., go < gi (see Table Such a spe c trum is 
similar to the theoretical predictions bv' Ptuskin et al.l (I20I0I) . 

Finally, in Calculation Si (discussed in Section IZZt all pa- 
rameters are the same as in Calculation R, except for the diffu- 
sion coefficient Dq and the Alfven speed vau- Dq is reduced 
by approximately 25%, whi ch makes the rat i o zh/D u equal 
to that in the calculations of iBlasi & Amatol ( l2012a l). vau is 
reduced accordingly to avoid a large bump in the proton spec- 
trum at low energies. In addition, the propagation calculations 
use a diff erent set of total inela stic cross sections (Equations 
[6]-[8] in lHorandel et al.ir2007h . Calculation S2 has the same 
parameters as Calculation Si, but the density of all gas com- 
ponents in the Galactic disk (i.e., H I, H II, and H2) is multi- 
plied by 2. The results of Calculation Si ^2 are compared with 
Calculation R in Figurel?] 

4. RESULTS 

The results of Calculations R, S, P, I, L, and H, as specified 
in Tables [T] and 121 are summarized in Figures l4l through [T4l 
Figures m and |5] show the proton and He spectra and their ra- 
tio, and Figure |6] shows the CR electron spectrum. Figur e [T] 
illustrating Hypothesis S, is discussed in detail in Section l272l 
and in the figure caption. 

Since the origin of the difference between the slopes o f the 
proton and He spectra was discussed in detail in Section IZ21 
we do not mention that topic in this section, instead concen- 
trating on scenarios explaining the spectral break. Calcula- 
tions P, I, L and H were designed to reproduce the observed 
proton. He and electron spectra and, therefore, cannot be used 
to constrain any of these scenarios. However, their predictions 
for CR anisotropy and the production of secondary species 
(B/C ratio, p flux, pip ratio, e+ flux and e+/ (e+ + e^) ratio) 
differ. These predictions are shown in Figures [SUlJ] Predic- 
tions for the diffuse Galactic 7-ray emission at intermediate 
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latitudes (10° < |6| < 20°) are compared with the data col- 
lected by the Fermi-LAT in Figure [14] 

4. 1 . Proton and He Spectra 

Proton and He spectra calculated for the different scenar- 
ios and their ratio are plotted in Figures |4] and |5] The bins in 
rigidity for protons are different from the bins for He in all 
experiments. Because of this, the experimental data points of 
PAMELA, ATlC-2, and CREAM shown in Figure|5]were ob- 
tained by interpolating the proton and He spectra, along with 
their errors, and by calculating the pfHe ratio on a grid. For 
simplicity, solar modulation for all spectra is taken into ac- 
count using the force-field approximation (G lee son & AxfordI 
[1968) with a modulation potential $ = 450 MV. 

While the reference case. Calculation R, provides satis- 
factory agreement with pre-PAMELA data by construction, 
it naturally misses all newly discovered features: an overall 
harder He spectrum, the spectral break at pbr, along with the 
dip just below pbr- The difference between the spectrum of 
He and protons at all energies was phenomenologically in- 
cluded in all other calculations except Calculation R, which is 
reflected in the considerably better agreement with the p/He 
ratio data in Figure [5] 

Calculation P. A break in the rigidity dependence of the 
diffusion coefficient leads to a corresponding break in the CR 
spectrum at Earth. In order to match the data, we assumed the 
change in the value of 6 (see Equation (O) from 6i = 0.30 to 
62 = 0.15 at pi ss pbi - For the PD model, the index changed 
at pi from 61 = 0.60 to 62 = 0.37. The difference 61 — 62 is 
greater in the PD than in the D-R model, because in the latter 
the reacceleration process additionally softens the spectrum of 
protons below po. Two corresponding physical quantities can 
be derived from these values. Assuming that the change in S 
is caused by a difference between the properties of interstellar 
MHD turbulence on scales smaller and larger than a certain 
length scale Abr, we can estimate this length to be of order 
of the gyroradius of 300 GV particles. The gyroradius of a 
particle of rigidity p in magnetic field B is 

For a characteristic interstellar magnetic field of order of a few 
pG, this implies a change in turbulence properties on length 
scales of the order of Abr ~ 10 AU. If the quasilinear the- 
ory of turbulent particle diffusion applies to CR transport in 
the ISM, the value 62 — 0.15 corresponds to turbulence spec- 
tral index, a — —2 + 62 = —1.85, which is harder than a 
Kolmogorov spectrum, a ~ —5/3. Note that the direction in 
which the index a changes across the transition wavenumber 
k = Ay^^ is opposite to the transition of the turbulent cas- 
cade from the inertial to dissipative regime. In our case, the 
turbulence spectrum must harden, rather than soften, above 

Calculation I assumes a change of the power-law index of 
CR injection spectrum at pi^ji — 300 GV, which produces a 
break at pbi ~ Pinji in the CR spectrum at Earth. The dip is 
not produced in this calculation. 

Calculation L (dashed orange lines in Figures [4] and O 
agrees with the p/He ratio and spectral break, and also repro- 
duces the dip just below pbi - This is possible because of the 
combination of the hard spectrum from Galactic sources that 
matches the data for p > pbr (solid orange lines in Figure |4]i 



with a local low-energy source having a sharp turnover just 
below pbr (dotted orange lines in Figure |4|i. 

In Calculation H, the spectral break at pbr is produced 
by the local source beginning to dominate the CR spectrum 
above pbi - We assume that the Galactic source has the same 
power-law index for p > pbi as the low-energy CR spectrum. 

The observed continuity of the p/He ratio and its slope at 
Pbr within statistical and systematic uncertainties is very im- 
portant. In Scenario P, this property of the p/He ratio comes 
about naturally. Indeed, if the injection spectrum is continu- 
ous, then protons and He nuclei experience the change in dif- 
fusion coefficient in the same way, and the p/He ratio is unaf- 
fected. However, matching this observation in the framework 
of a composite source spectrum (Scenario L, Scenario H, or 
Scenario I (&)) requires an additional assumption of the H 
to He ratio to be the same at the sources producing the low- 
energy and high-energy particles. 

For the analysis of all calculations discussed above, the dip 
in the spectrum, if it is significant, may lead to important im- 
plications. One possible explanation for the dip may be pro- 
vided in the framework of Scenario L and Scenario I (6). 
It can naturally appear if the spectrum of the low-energy 
CR sources (local, as in Scenario L, or Galactic, as in Sce- 
nario I (b)) sharply turns over just below pbr, rather than con- 
tinuing as a power law up to the knee in the CR spectrum. In- 
deed, it is trivial to prove that for any two power-law spectra, 
their sum always hardens with energy. Thus, for the softening 
below pbi to occur, the low-energy source spectrum may not 
be a pure power law; the sharpness of the dip suggests that it 
must steeply turn over below pbi , where the dip occurs. The 
dip may also be explained in the framework of Scenario P, if 
a corresponding dip in the spectrum of MHD turbulence re- 
sponsible for CR confinement in the Galaxy is assumed. It is 
not possible to explain the dip with Scenario H because the 
low-energy source is assumed to have a power-law spectral 
shape extending all the way to the knee. 

4.2. Electrons 

The CR electron spectrum in this problem is connected to 
the proton spectrum because (1) electrons propagate in the 
Galaxy in the same magnetic fields as nuclei, and (2) some, 
if not all, CR electrons are produced by the same sources as 
nuclei. 

As Figure [6] shows, the propagated electron spectrum in 
Calculation R does not fi t the observatio ns of the Ferw i-LAT 
(lAckermann etalJ l2010^ or PAMELA (lAdrianillMTI) . In- 
deed, the observed spectrum appears convex (i.e., harden- 
ing with increasing particle energy) from wlO to 10"^ GeV, 
whereas the calculated spectrum in this energy range is con- 
cave. The concavity is caused by energy losses on ionization 
at low energies, and synchrotron losses at high energies. 

In Calculation P, despite a break in the diffusion coeffi- 
cient, the electron spectrum at Earth does not fit the high- 
energy data and is not convex (Figure |6]l. This is because 
synchrotron energy losses above 100 GeV oppose the effect 
of the diffusion coefficient break, and cause the spectral soft- 
ening. 

In the injection effects Scenario, represented by Calcula- 
tion I, it is possible to modify the source spectrum of electrons 
in order to fit the data. Indeed, since the nucleon injection 
spectrum has a break at pbr, it is natural to assume that the 
electron injection spectrum may have a similar feature. This 
argument holds both in the source with a spectral break in- 
terpretation (Scenario I (a)) and the composite spectrum in- 
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terpretation (Scenario I (&)) of Calculation I. Moreover, the 
energy and magnitude of the break do not have to be the same 
for electrons and protons because the electron to proton ra- 
tio may vary with energy and with source type. As Figure |6] 
shows, one can achieve agreement with the data by assuming 
that the electron source spectrum has a break (hardening) at 
i? = 70 GV (61 GV in the PD model), with the index chang- 
ing from -2.70 to -2.33 (or -2.67 to -2.24 for PD). Note 
that one cannot justify such a break in any other scenario, un- 
less CR electrons and nuclei are assumed to be produced by 
different sources. 

In Calculation L, we assumed that above the break (> 
100 GV for electrons) the particles are produced by a hard 
Galactic source, and that below the break the electron flux is 
dominated by an unknown /ow-energy local source. Likewise, 
in Calculation H, an unknown high-energy source of CR elec- 
trons was assumed. In both cases, the flux of the local source 
was calculated as the difference between the observed and the 
calculated Galactic source fluxes. Therefore, the total spec- 
trum agrees with the data at all energies. Note that only the 
Galactic source flux was used in the calculation of secondary 
lepton production and 7-ray emission. 

4.3. Anisotropy 

For all scenarios, we calculated the anisotropy of the high- 
energy CR flux at the location of the Sun due to diffusive es- 
cape of CRs from the Galaxy. The results are presented in Fig- 
ure [8l along wit h data. References to individual experiments 
may b e found in lPtuskin et alj (l2006al) : see also IStrong et alj 
(120071) for a color version of the plot. The anisotropy, dom- 
inated by the radial component, is highly sensitive to the 
choice of the diffusion coefficient and the spatial distribution 
of CR sources. 

Our calculation ignores the eff ect of nearby CR 
sources, which may be significant dPtuskin et alj l2006at 
iBlasi & Amatol l2012bl) . but is not very well defined and 
depends on the assumed distances to the sources and their 
ages. However, if the diffusive component of the anisotropy 
dominates in a certain energy range, two conclusions may 
be drawn from this plot. The first one is that Scenario P 
can be distinguished from the others with improved CR 
proton anisotropy data. The second point illustrated by our 
calculation is that the diffusion regime of Calculation P with 
D oc p^^ and ^2 = 0.15 agrees with the available data better 
than §2 = 0.30. The PD model, due to a harder dependence 
of the diffusion coefficient on rigidity, predicts a higher 
degree of anisotropy, which disagrees with the data more than 
the other calculations. 

In the case of Calculation H, the plotted lines correspond to 
only the Galactic source, while the direction and magnitude of 
the local source flux anisotropy above pbr are unknown. De- 
pending on the location and proximity of the local source, the 
direction of the overall CR drift can be changed substantially, 
because the local source flux above 10"^ GV is comparable to 
the Galactic flux (see the bottom plots in Figure H]). Since 
the estimate of the distance to the local source is beyond the 
scope of this paper, we do not provide quantitative predictions 
of anisotropy in this case. 

4.4. Boron to Carbon Ratio 

The B/C ratio for all scenarios discussed in the paper is 
shown in Figure |9l Predictions of Calculation R and Calcu- 
lation I coincide at all energies, while Calculation P predicts 



a larger B/C ratio for p > pbi , which is a consequence of the 
smaller diffusion coefficient in Calculation P. In the case of 
Calculation L and Calculation H, the results include the effect 
of the local CR source. CR boron is produced by fragmenta- 
tion of heavier elements and decay of ^"^Be. If the local source 
is very nearby, its flux should contain no boron. However, the 
abundance of (primary) carbon should be close to the inter- 
stellar value, which results in a lower B/C ratio in the net flux 
than without the local source. 

To find the B/C ratio for Calculation H and Calculation L, 
we assume that the local source produces no boron and as- 
sume that the flux of local source carbon is proportional to 
that of He, with the same carbon and He abundance as in 
the Galactic source. Calculation H predicts a lower B/C for 
p > pbi because the local high-energy source supplies pri- 
mary carbon, but not secondary boron, at these energies. Cal- 
culation L could not be tuned to reproduce the B/C ratio in the 
D-R model (see Section [33] ) However, in the PD model, the 
only calculation fitting the B/C ratio below 1 GeV nucleon^^ 
is Calculation L due to the reduction of the B/C ratio by the 
local source. 

Experimental data at low energie s (below 
1 GeV nucleon^^) were collected by ACE (iDavis et"ar 
2000), and for hig h energies by HE AO-3 (jEngelmannetaL 
199Q), CREAM ("Ahnetal.l I200 I, ATIC-2 (.Panov et al. 
20081) . and TRACER ( Obermeier etal] 1201 Ih. The uncer- 
tainties in the data are still too large to rule out any of the 
scenarios considered in this paper, but data collected by future 
experiments, such as AMS-2, may be more constraining. 

4.5. Antiproton Flux and p/p Ratio 

The antiproton flux, another probe of CR propagation, is 
plotted in Figure [TOl and p/p ratio in Figure [TTl together 
with the PAMELA dat a from lAdriani et all ( 120101) and the 
BESS-Polar II data from' Abe et al.l(l20Ilh . Calculations R, P, 
I, and H are in good agreement with data below 100 GeV. Dif- 
ferences between all calculations are apparent above ~1 TeV, 
but no data are currently available. 

Calculation L, the only case of the PD model where the 
low-energy B/C data are reproduced, predicts a factor of '--^ 2 
excess of p below 100 GeV, due to a larger particle confine- 
ment time. In Calculations L and H, the local source was 
assumed to be completely devoid of primary or secondary an- 
tiprotons. More accurate data covering a larger energy range 
may help eliminate some of the scenarios, and the AMS-2 
mission may provide these data. 

4.6. Positrons 

Figure[T2]shows the calculated positron flux, and Figure [T3] 
— the positron fraction (e+/(e~ + e+)). Our model does not 
include a source of primary CR positrons; the e+ particles 
in all our calculations are produced in inelastic collisions of 
other CR species. 

The positron flux measured by the Fermi-LAil 
dAckermann et al.l 120 12h is significantly greater than the 
prediction of all our calculations (Figure [T2] |. The only 
scenario that could explain the discrepancy is Scenario H in 
which the local high-energy CR source produces primary 
positrons. 

The positron fraction (Figure [T3T l provides additional 
evidence that the calculations predict insufficient flux of 
positrons at high energies. However, below a few GeV, the 
positron fraction in the D-R model is overpredicted. 
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4.7. Diffuse "f-Ray Emission 

Predictions of the 7-ray emission at intermediate Galactic 
latitudes (10° < |5| < 2 0°) are plott e d in F igure [141 along 
with the data reported by lAbdo et aTl (|201(^ see the o nline 
supplementary material). Following lAbdo et al.l ( 120101) . the 
flux of the inverse Compton (IC) component was increased by 
a factor of two to obtain a good fit to the data. Our calculations 
include 7-ray emission produced by hadronic and leptonic 
components of CRs (i.e., 7r°-decay, IC, and bremsstrahlung 
channels), as well as point sources, and the isotropic extra- 
galactic emission. The relative differences in the total 7-ray 
flux between the considered scenarios are quite small and are 
considerably smaller than if only the 7r°-decay channel is con- 
sidered (as in Donato & Serpico 2013]). 

Predictions of all calculations, except Calculation L, agree 
with the published Fenni-LAJL data, within the uncertainty 
band. Calculation L predicts a slightly lower 7-ray emission 
below 10 GeV. Note that even the reference Calculation R, 
which does not agree with the PAMELA data, satisfactorily 
reproduces the 7-ray data. 

For all scenarios we calculated the 7-ray spectrum up to 
1 TeV, but the Ferm/-LAT team has not published on the data 
above 100 GeV so far. At these energies, the softer spec- 
trum of protons (above phi) from Galactic sources in Calcu- 
lation H produces less pions resulting in a smaller flux of pio- 
nic 7-rays compared to other scenarios. However, the contri- 
butions of comparable IC component and isotropic emission 
in the range 100 GeV - 1 TeV are not affected by the pro- 
ton spectrum. Therefore, even at these energies the difference 
between the total 7-ray emission in Calculation H and other 
calculations is significantly smaller than the difference in the 
7r"-decay channel alone. Unsurprisingly, Calculation P can- 
not be distinguished from Calculation I using the 7-ray data 
alone, because calculations for both scenarios result in nearly 
the same spectrum of CR protons, even though it is achieved 
via different mechanisms. 

As the Fermi mission continues, the statistical uncertainty 
will be reduced as data accumulates, and systematic errors are 
likely to be brought down by improved data analysis. It should 
be noted, however, that the analysis of the diffuse 7-ray emis- 
sion is complicated by many factors, including the uncertainty 
in the spatial distribution for the CR sources and the loosely 
constrained spectrum of CR electrons over the Galaxy respon- 
sible for IC emission that dominates high-energy 7-rays. 

5. SUMMARY 

We have presented scenarios reproducing the spectral fea- 
tures in CR proton and He spectra (the p/He ratio dependence 
on energy, the dip, and spectral break) observed by ATIC- 
2, CREAM, and PAMELA. For each scenario, we performed 
CR propagation calculations in the framework of the D-R 
model (except Scenario L), using the GALPROP code. Dif- 
ferences between scenarios are reflected in the CR anisotropy 
and fluxes of secondary CR species: the B/C ratio at high en- 
ergies, the antiproton flux and antiproton to proton ratio, as 
well as the diffuse Galactic 7-ray emission. We find the fol- 
lowing: 

(1) He spallation (Hypothesis S) may be partially responsi- 
ble for making the spectrum of He and heavier nuclei harder 
than protons. However, a significantly increased grammage 
traversed by CRs in the Galaxy is required to explain the p/He 
observations with spallation alone. This makes it problematic 
to match stable secondary CR isotope observations (B/C and 



antiprotons). 

(2) Electron spectrum can be reproduced in Scenario I (a) 
(break in the injection spectrum) or Scenario I (b) (compos- 
ite Galactic source) only if the break in the electron injection 
spectrum is stronger, and occurs at a lower rigidity, than in the 
proton spectrum. A break in the diffusion coefficient (Sce- 
nario P) cannot simultaneously explain the concavity of the 
observed proton and electron spectra. 

(3) Experimental uncertainty in the data on high-energy 
B/C ratio does not allow us to rigorously reject any of the 
scenarios for the origin of the spectral break. However, more 
accurate measurements of high-energy B/C, expected from 
planned CR experiments, may be used for model rejection. 

(4) In the D-R model, low-energy B/C data are consis- 
tent with any scenario except Scenario L (low-energy local 
source). In the PD model, low-energy B/C data require that 
the diffusion coefficient decreases with increasing energy; 
however, it is also possible to achieve agreement with these 
data in Scenario L for constant diffusion coefficient. 

(5) Antiproton flux and p/p ratio seem to disfavor the low- 
energy local source hypothesis (Scenario L). Measurements 
of p and/or pip above 1 TeV may help to differentiate between 
the other scenarios. 

(6) Radial component of the diffusive anisotropy of CR flux 
is too high in all scenarios, but the discrepancy is larger in Sce- 
nario L, while Scenario P predicts the lowest anisotropy. The 
PD model predicts a higher anisotropy than the D-R model. 
Local sources may significantly affect the CR anisotropy, and 
therefore our simple analysis applies only to energy range un- 
affected by local sources. 

(7) Data on the positron flux and positron fraction are incon- 
sistent with any of the scenarios, if all observed CR positrons 
are secondary. However, if some of the detected positrons are 
produced in sources, then only Scenario H (local high-energy 
source) can account for the observed positron excess at high 
energies. 

(8) Finally, the 7-ray data are in agreement, within the un- 
certainty range, with all scenarios, including Scenario R, even 
though the reference scenario does not agree with the new 
measurements for the CR proton and He spectra. Scenario L 
slightly underpredicts the 7-ray flux below a few GeV. 

Most specific physical models explaining the p/He ratio, 
spectral break and the dip fall into one of the scenarios stud- 
ied in this paper, or their combination. Data from experiments 
such as the ferm/-LAT and AMS-2 can be used to distinguish 
between some of these scenarios. 
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Table 2 

Summary of Model Parameters, Plain Diffusion (PD) Model 
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Figure 1. (color in online version) Galactic CR source injection spectrum for Calculations R, S, P, I, L, and H in arbitrary units. Left: diffusive-reacceleration 
model, right: plain diffusion model. The normalization for the injection spectrum was chosen to match local measurements of proton and He spectra. For all 
calculations, the lines represent the Galactic CR source injection spectrum. "Local" sources, present in Calculation L and Calculation H, are not shown here. 
The "local" source fluxes at Earth in Calculation L and Calculation H were obtained as the difference between the observed and propagated Galactic fluxes. 
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Figure 2. (color in online version) Galactic CR source electron injection spectrum for Calculations R, S, P, I, L, and H in arbitrary units. Left: diffusive- 
reacceleration model, right: plain diffusion model. The normalization for the injection spectrum was chosen to match local measurements of proton and He 
spectra. For all calculations, the lines represent the Galactic CR source injection spectrum. "Local" sources, present in Calculation L and Calculation H, are not 
shown here. The "local" source fluxes at Earth in Calculation L and Calculation H were obtained as the difference between the observed and propagated Galactic 
fluxes. 




Figures, (color in online version) Diffusion coefficient of CRs in the Galaxy (see Equation (5)). Left: diffusive-reacceleration model, right: plain diffusion 
model. The values of the D coincide for Calculations R, L L and H. For Calculation S, the value of D is slightly smaller at all energies. For Calculation P, a 
break in the diffusion coefficient is assumed, changing 81/82 = 0.30/0.15 at pi = 300 GV. 
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Figure 4. (color in online version) Propagated CR pr oton and He spectra: data of PAMELA <Adriani et ZII201 ID , ATIC-2 jWefel et al]|2008t IPanov et al] 
I20Q91) . and CREAM jAtin et aIll201(]trYoon et al.ll2011D together with calculation results. Left: diffusive-reacceleration model, right: plain diffusion model. For 
Calculations L and H, solid lines show the net CR flux comprised of the Galactic source contribution (dashed lines) and a "local" source contribution (dotted 
lines). Sola r modulation in all spectra is taken into account using force-field approximation with a modulation potential $ = 450 MV. See the discussion in 
SectionRTT] 
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Figure 5. (color in onlin e version) Propagated C R proton to He flux ratio: data of PAMELA I Adriani et al. 2011), ATIC-2 i Wefel et al. 2008; Panov et al. 2009?), 
and CREAM I Ahn et al. 2010t rVbon et al.ll201 II) , together with calculation results. Left: diffusive-reacceleration model, right: plain diffusion model. The p/He 
points for experimental data were o btain ed by interpolating the measured fluxes of protons and He along with respective eiTors and calculating the p/He ratio on 
a grid. See the discussion in Section l4~l1 
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Figure 6. ( color in online version) P ropagated CR ele ctron spectra: m odels and data. Left: diffusive-reaccelerat ion model, right: plain diffusion model. The 
data are from lAckermann et afl 120 101) (Fermi-LPJ:) and lAdrianil 1201 II) (PAMELA). See the discussion in Section|42] 
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Figure 7. (color in online version) CR proton to He, B/C ratio and antiproton flux (sources of data as in Figuresl5]|9l llOt . together with calculation results. Left: 
diffusive-reacceleration model, right: plain diffusion model. Calculation Si and Calculation S2 use the same spallation cross sections and / D ratio as those 
of Blasi & Amato (2012a). In Calculation S2, the gas density is increased by a factor of 2 with respect to Calculation Si . Results of Calculation Si are consistent 
with the B/C ratio data, but slightly overpredict p measurements. However, note that the effect of spallation is not sufficient to make the slope of the p/He ratio 
agree with the PAMELA data. In Calculation S2 , spallation is stronger due to increased grammage, and the p/He ratio is reproduced better, but the B/C ratio and 
p flux are overpredicted. 



Testing the origin of HECRs 



19 




lO' lO"* 10" 10** lO' lO" 10^ lO'* 

E, GeV E, GeV 

(a) (b) 

Figure 8. (color in online version) CR flux anisotropy: data from'Ptuskin et al. (2006a'), together with calculation results. The anisotropy was calculated as the 
ratio of diffusive flux in the radial direction to the isotropic flux at the corresponding energy. Left: diffusive-reacceleration model, right: plain diffusion model. 
Calculations R, S, I, L and H all predict the same value of anisotro py, a nd their respective lines overlap. The result of Calculation P is different due to a different 
form of the CR diffusion coefficient. See the discussion in Section l43] 
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Figure 9. (color in online vers ion) CR bo ron-to-carbon flux ratio. Left: diffusive-reacceler ation model, right: plain diffusion m odel. Data: IDavis et al] 120001) 
(ACE). Engel mann et al.'(T9^ (HEA0-3).'Ahn et aUnwB) (CREAM).! Panov et aUHW^ (ATIC-2) and lObeniieier et aI.I<20TTI ) (TRACER). For Calculation L 
and Calculation H, dashed lines show the ra tio o f just the Galactic source, while solid lines show the B/C ratio including the contribution of the "local" source 
component. Additional discussion in Section l44l 
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Figure 10. (color in online version) CR antiprotons; data fiom lAdriani et al] j2010l) (PAMEL A) an d lAbe et al.1 )201 II) (BESS) together with calculation results. 
Left: diffusive-reacceleration model, right: plain dilfusion model. See the discussion in Section l431 
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Figure 11. (color in online version) CR antiproton to proton ratio: data from A driani et all j201QD (PAMELA) together with calculation results. Lefi: diffusive- 
reacceleration model, right: plain diffusion model. See the discussion in Section [431 
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Figure 12. (color in o nline version) Positron flux: models and d ata. Left: diffusive-reacceleration model, right: plain diffusion model. The data are from 
lAckermann et alj j2012l) (Fermi-hAI). See the discussion in Section l431 




Figure 13. ( color in online version) Positron fraction: models and data. Left: diffusive-reacceleration model, right: plain diffusion model. The data are from 
lAdriani et al] (2009.) (PAMELA) and Arkermann et al. (2012) (feraii'-LAT). See the discussion in Sectionl431 
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Figure 14. (color in online version) Diffuse 7-ray emission from intermediate Galactic latitudes: data from I Abdo et al J 420100 (Feraii'-LAT) together with 
calculation results. Left: diffusive-reacceleration model, right: plain diffusio n model. The dat a are the diffuse 7-ray intensities averaged over all Galactic 
longitudes and intermediate Galactic latitudes 10° < |fe| < 20°, as reported by Abdo et alj RCTTO (available in the online supplementary material to the article). 
See the discussion in Section l477l 



